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An exact formula is given for the probability that there exists a spanning cluster between opposite 
boundaries of an annulus, in the scaling limit of critical percolation. The entire distribution function 
for the number of distinct spanning clusters is also derived. These results are found using Coulomb 
gas methods. Their forms are compared with the expectations of conformal field theory. 



, Since Langlands et al. [1] conjectured on the basis of numerical evidence that crossing probabihties between two 
' non-overlapping segments of the boundary of a simply connected region should be conformally invariant, there has 
, been intense interest in the scaling limit of two-dimensional percolation [2]. In [3] it was shown that this invariance 
^' was implicit in ideas of conformal field theory, which in addition yielded an explicit formula. Further exact formulas 
^ ^ were conjectured, by Watts [4] for the probability of a simultaneous left-right and up-down crossing, and by Pinson 
C/^ ' [5] for various crossing probabilities on the torus. The latter work used so-called Coulomb gas methods [6], which had 
, been developed for more general two-dimensional critical systems, in parallel with those of conformal field theory. In 
' [7] results were conjectured for the asymptotic behaviour of the probabilities that at least Nc distinct clusters cross 
either a rectangle or an annulus, using earlier conjectures of Saleur and Duplantier [8]. In [9], among other results, a 
\l ' prediction was given for the mean number of crossing clusters in the opposite limit, when this number is large. 

Meanwhile, starting from another approach, Schramm [10] conjectured that the scaling limit of percolation hulls is 
generated by stochastic Loewner evolution (SLEg). From this follow many results [11], including the original crossing 
formula and the exponents in [7]. Finally Smirnov [12] proved the original crossing formula for site percolation on the 
00 triangular lattice, and hence the validity of the SLEg approach to percolation [13]. 

In this paper we refine the results of [7] for the annulus, presenting results for a general value of the modulus. 
Consider a critical percolation problem in a non-simply connected region of the plane with the topology of an annulus. 
The boundaries are assumed to be suitably smooth. The interior of this region may be conformally mapped into the 
interior of a circular annulus Ri < \z\ < i?2, with modulus q = {R1/R2), or into the rectangle {Q < x < £,0 < y < L) 
with the edges at a; = and x — £ identified, and q — e"^^^/^. A crossing (or spanning) cluster is one which contains a 
' path connecting the opposite boundaries. Let P{Nc) be the probability that there are exactly Nc non-overlapping such 
"j^ I clusters. When A^c = 1 it is possible for the cluster also to wrap around the x-cycle on the annulus. By convention, 
' we do not count such clusters as spanning. 
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Results. 



^ The crossing probability is 



X! -P(^c) - %/3 ('~i2r-2+2.r ^12-r2 + 10r+2\ 



Furthermore we have an explicit expression for P{Nc) for Nc > 1: 



PW - Iin:^ n (1 - fn-'T.MN.) g^-A (2) 



where 



sj \ 2r 

These results may be transformed into other expressions using various theta-function identities. For example, in terms 
of the conjugate modulus q = e~^^/^, we find 
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jVc=i Z^rez l^y y ; 

Note that for L/£ > i only a few terms need be kept in (1,2) for great accuracy, while for L/£ < ^ the same is true 
of (4). 

Both the numerator and denominator of (1) and (4) are specialisations of Jacobi theta functions, and hence may 
be written as infinite products. In terms of Dedekind's eta function r]{T) = 9^^^'* nJ^Li(l ~ 9")> with q = e"^^^^ , we 
find 

" %(3r)r?(2r)2 " ^(-l/3r) r?(-l/2r)2 " 



Coulomb gas method. 

Although wc shall later argue that these results arc indeed conformally invariant, it is simpler first to set the problem 
up in the periodic rectangle defined above. Consider a portion of a regular triangular lattice covering the rectangle, 
oriented as shown in Fig. 1, so that, if the lattice spacing is a, there are 2{£/a) + 1 columns and (2/-\/3)(L/a) + 1 
rows of the lattice. Impose periodic boundary conditions in the x-direction, so that the rightmost column is identified 
with the leftmost one. Consider a critical site percolation problem on this lattice, in which sites are independently 
coloured red or blue with equal probability. A (blue) cluster is a set of blue sites in which every site is connected 
to every other by a path which traverses only blue sites. A spanning c;lust(a' is one which contains at least one site 
on the edge y = and one site on the edge y = L. For a particular assignment of colours, let be the number of 
distinct non-overlapping spanning clusters. We are interested in the distribution P{Nc) of the random variable Nc in 
the continuum limit a ^ 0, for fixed L and Scale invariancc implies that it should depend only on C./L. 

Instead of considering the clusters, we may equivalently consider the configuration of hulls which separate them. 
A hull is a path on the dual lattice (in this case a honeycomb lattice) which separates blue sites from red sites. In 
our case, hulls can either form closed paths or be open, each end terminating at an edge. Open hulls which have 
ends terminating on different edges are called spanning hulls. The number of such spanning hulls is ric = 2Nc. An 
allowable configuration of hulls is one in which each dual site is connected to either or 2 neighbouring dual sites, 
except for the edge sites, which may be connected to cither or 1 neighbouring site. In addition, the number of 
spanning hulls must be even. The correct weights are achieved by weighting all allowable hull configurations equally. 

A related model is the 0(1), or Ising, model on the dual lattice, at zero temperature. In this model Ising spins 
s(r) = ±1 reside at each site r of the honeycomb lattice, and the partition function is 

^0(i)=n E n Kl+Mr)s(r')) (6) 

r s(r)=±l (r,r') 

where the latter product is over all nearest neighbour pairs (r, r'). The 'high-temperature' expansion in powers of 
i, afterwards setting i = 1, reproduces exactly the hull configurations of the percolation model in the case when the 
number of spanning hulls, denoted by n,;, is even and equal to 2Nc, but in the 0(1) model Uc may also be odd. 
Evidently when i = 1, ^o(i) = 2. Denoting by p{nc) the probability that there are exactly ric spanning hulls, we 
have P{Nc) — 2p{2nc). The first factor of 2 arises because to each allowable configuration of hulls there correspond 
two assignments of colours. We can construct the generating function by weighting each spanning hull by a factor u. 
Denoting the corresponding partition function by Z{u) we therefore have 

oo 

^ p{n,)u-^ = \Z{u) (7) 

nc=0 

oo 

^ P{N,)u^^^ = \{Z{u)^ Z{-u)) (8) 

Each hull may assigned a random orientation, so that to each configuration of H hulls correspond 2^ configurations 
of oriented hulls. The weights for each orientation should be chosen so that they sum to unity (resp. u) for each 
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(spanning) hull. For closed hulls, it is conventional [6] to assign a 'weight' e^*'^ to each dual site at which an oriented 
hull turns through an angle of ±7r/3, where x = tt/IS is chosen so that the total weight for a closed hull, on summing 
over its orientations, is e^^^ +e~^'^^ = 1. However, this does not correctly account for closed hulls which wrap around 
the x-cycle of the annulus, which would have weight 1 + 1 = 2 according to this scheme. Such configurations can only 
occur when Uc = 0. Thus, for the time being, we assume that ric > 1. The case Uc = may be inferred afterwards 
using the overall normalisation of the partition function. 

For oriented hulls which terminate at an edge, let us assign the same weights as above for internal turnings, and in 
addition weights a or /3 to their extreme segments, as shown in Fig. 1. By choosing 



^ . e'^^r^ and /? = e'^''^ , (9) 

cos 3x7 V cos 3x7 

hulls which begin and end on the same edge are counted with a weight (cos6x/ cos3x)(e3*'' +e~3»x^ = \^ as required, 
while spanning hulls carry a weight (cos6x/ cos3x)(e^''^ + e~3»x )^ go that we should identify 

■u = cos 3x'/ cos 3x • (10) 

The factor (cos6x/ cos3x)^^^ coming from (9) is raised to a power E which is the total number of ends of open 
hulls, whether they be spanning or not. An open end occurs every time the neighbouring sites of the triangular 
lattice are of opposite colours. Since these arc independently distributed, i;^ is a sum of 0{2£/a) independent"'^ 
random variables, each taking the values or 1 with equal probability. In the continuum limit a/^ — » 0, therefore, 
(cos6x/cos3x)^^^ ~ (cos6x/cos3x)^^^", with probability one. These contribute to the non-universal edge free 
energy, but not to the universal dependence on i/L. 

Let Z{3x, 3x') be the partition function of the loop gas with the above phase factors but ignoring the factors of 
(cos6x/cos3x)^/^. Then 

OO 

^ p(ne)u"^ = Ci (Z(7r/6, xO - Z{7r/6, it/2)) (11) 

nc = l 

where Ci is a non-universal number and the second term, with cos3x' = 0, subtracts out the contribution with ric = 
which is incorrectly counted by the above scheme. 

The configurations of the oriented loop gas are in 1-1 correspondence with those of a height model on the original 
triangular lattice. These heights h{r) are conventionally chosen to be in ttZ, and are defined by the conditions that 
/i = at some fixed site, say (0,0), and that it increases (decreases) by ±7r each time an oriented hull segment is 
crossed. On the annulus, however, we must also allow for possible jumps A/i > [7r| in h across some path which spans 
the annulus, say along x = — \a. The factors e^^^x /2 then accumulate to e^^^ Qg^]x edge. 

So far, everything is finite and exact. In the conventional Coulomb gas method [6], one now assumes that in 
the continuum limit [aj l^aj L) — s- we may replace /i(r) by a rc^al- valued field, with a gaussian measure cx exp ( — 
(Sf/47r) / {dh)'^dxdy) . For the models we axe considering, g is fixed to be 1 — (Sx/tt) = |. We shall assume that the 
same is true on the annulus, except that we must allow for a possible discontinuity around the a;-cycle. Thus we write 

h{x,y) = {pTT/e)x + h{x,y) (12) 

where p £ Z and h{x + i,y) = h{x, y). Substituting in this decomposition, 

where 2 (x J X)he~^^^'^^-' ■f^^'^'' '^'^'^^ is the universal part of the partition function of a free field on the annulus, with 
Neumann boundary conditions, and with the constraint that h{Q) = 0, which removes the zero mode. The factor C2 
is non-universal, and reflects the contribution of the short-distance degrees of freedom which are integrated out in the 
coarse-graining assumed in adopting the gaussian measure. It is expected to depend exponentially on the total area 
(tLja^) and the perimeter {2i/a), but is not expected to have nontrivial dependence on the modulus l/L. 



^Almost independent, since the sum along each edge must be even. 
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The c = 1 partition function Z is well-known [14]. Writing it as Tre~^^^, where Hh is the quantum hamiltonian 
for a free field on circle of perimeter L, it is H^i X]5/=o e~^-^"'" where -E„,]v = (A'' + 5)(n7r/L). The leading term as 
IjL ^ CO comes from N = and is proportional to HJ^i e"''^^/^^)". This must be regularised, however. Apart from 
a cut-off dependent term which can be absorbed into C2, it gives e"^'^^/^^^^^""'^-' = ■ where q = er'"^l^ . The terms 
with iV > 1 give Jli^iCl ^ 9")^^- The zero-mode h = constant is suppressed in the functional integral over h since 
we set /i.(0) = 0. However going from this constraint to one on the n = mode introduces a jacobian proportional to 
(i/£)V2 [14]. Finally we have 

00 

n=l 

Eq. (13) may now be transformed using the Poisson sum formula: 



rez * 



rt=l r6Z 

Note that the {L/iY^"^ factors cancel. 

Setting now 5 = § and subtracting the contributions with 3%' = f and 3%' = f , we arrive, after some algebra, at 
the result for the 0(1) model 

Because our height model phase assignments incorrectly count loops which wrap around the 3;-cyclc, we cannot 
directly compute the contribution with Uc = and therefore cannot fix C5 by demanding that '^^^^Qpijic) = 1- 
However, since =qP{^c) = 2X^jv pC^^c) = 1, it follows that Yl'^^^oPi'^n + 1) = |, and we can compute this in 
terms of Z{u) - Z{-u) oc Z{-k/Q,'k/&)- Z{-K/Q,b'Kl&). This gives 

S /'g6r^4-r _ ^6r^-|-5r-l- A 

E P^^o) = \C, "^^1^^ = \C, (18) 

odd lk=iU 9 J 

where the last equality follows from Euler's pentagonal number theorem [15]. We conclude that C5 = 1. (17) is then 
our main result for the probability of a crossing in the 0(1) model. 

For percolation, we need to compute =iP(2A^c) oc Z(7r/6, -k /&) + Z{-k /&, 5it/6) — 2Z(7r/6, 7r/2). This gives the 

main result (4). 

With the knowledge that C5 = 1, we may now transform these results back into series in g' = e"^'^^'/^ = g^'^'^. 
Using the identity r}{T) = (— ir)~ 2 jj(— 1/r), and the Poisson sum formula, we find 

00 

Z(./6, ^')='Ylil- ^ e^-'^ g^P'-^y^' (19) 

Thus for the crossing probability in percolation we have 

00 ^00 

P(7Ve) = n (1 - 9'")"' E (^°« f + COS 5|£ - 2 cos f ) (20) 

The expression in parentheses takes the value 3 if p = ±2 (mod 12), the value —3 if p = ±4 (mod 12), and vanishes 
otherwise. This leads to the first form (1) of our main result. The numerator in this expression may also be written 
as [16] 
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=gll?4(2T|6T) (21) 

ri6Z 

oc 

= n - - - ' (22) 

ra=l 

which, after a few more manipulations, gives (5). 

In order to find an expUcit formula for p{nc), we should solve (10) for e^'^^ in terms of u, which gives 

gSix' = e^-/2 (^^i _ (V3w/2)2)i/2 _ e-'-/2(V3w/2)) . (23) 

Substituting this into (19), expanding in powers of u, and identifying the coefficient of m^^", then leads to the result 
in (2,3). 

A further check on our results comes from differentiating (11) with respect to u at u = to obtain the mean number 
of crossing clusters. In the limit £ ;» L wc find E[Nc] (\/3/4)(^/i), in agreement with [9], and with a rigorous 
result of Smirnov [12] for the triangular lattice. 

If instead of the periodic rectangle we have a more general annular region, in order that spanning and non-spanning 
hulls be counted with their correct weights and 13^ in (9) must be modified by factors e^'^/®, where 9 is the (signed) 
angle which the tangent vector at the boundary makes with the x-axis. However, since the boundaries form simple 
closed curves, these extra factors accumulate to unity on each edge. The calculation then proceeds as before, yielding 
a conformally invariant result.^ 



Relation with conformal field theory. 

The crossing probability ^iP{nc) in the 0(1) model may be expressed as a difference — of partition 

functions in the n 1 limit of the 0{n) model, where Z_| denotes the partition function with the spins fixed in 

different directions on opposite edges, and with them fixed in the same direction. Evidently = 1, so that, 
from (17) 

= Ff35 — Tj zr{ (24) 

According to general boundary conformal field theory (BCFT) [18], any partition function like this should be ex- 
pressible in the form dhq^ , where h runs over all boundary scaling dimensions and dh is a degeneracy factor. For 
unitary conformal theories this must be a non-negative integer, but this need not be true here. From (24) we identify 
the smallest value of h to be i; this is identified in BCFT [18] as the scaling dimension of the 'boimdary condition 
changing operator' (/)+|_. This is consistent with the analogous result for percolation: see [3]. From (24) we see there 
is also an operator with h= 1. This we tentatively identify as introducing a hull which wraps around the a;-cycle in 
the 0(1) model, but does not touch either edge. This should carry an 0(n) index c which is not equal to either -I- or 
— , otherwise it could be absorbed at the edges. There arc n — 2 = — 1 possibilities for c, which accounts for the fact 
that this state occurs with degeneracy (—1) in (24). 

The powers (4p2 — 1)/12 in (2) are the well-known bulk multi-hull dimensions for percolation [8] . In accordance with 
general ideas of BCFT [18], (24) may be written as a linear combination of Virasoro characters Xh{<l) of irreducible 
representations of highest weight h, and, equivalently, as a combination of characters Xh{<l^), related by a modular 
transformation. This affords an interesting example of how BCFT works in a non-minimal theory. The details will 
be described elsewhere [19]. 



^Under a scale transformation r Ar, a partition function Z behaves in general as X'^^^^, where c is the conformal anomaly 
number, and, for a smooth boundary, x is the Euler number [17]. The latter vanishes for the annulus. If there are points 
on the boundary where it is not difFerentiable, however, there may be additional contributions [17]. In our case, these cancel 

between the c = 1 partition function Z and the Coulomb energy of the charges which accumulate at these singularities. This 
cancellation is connected with the fact that the overall conformal field theory has c = 0. 
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Comparison with numerical work. 



Shchur [20] has estimated -P(iVc) for A^c = 2, 3 and L = i. For this value o{ q = e it is sufficient to keep only 
the first term in (2), P{Nc) « 3^=-5g(4^c-i)/i2, which gives 

P(2) = 2.02 ... X 10"^ (exact), 2.0(4) x 10"^ (measured); (25) 
P(3) = 1.71 ... X 10"^ (exact), 1.4(5) x 10"^ (measured). (26) 

Our exact predictions fall within the (admittedly rather large) error bars. 



Summary. 



We have given explicit results for the probability that Nc critical percolation clusters cross an anmilus. From the 
point of view of conformal field theory, these results are different from the original crossing formula [3] in that they 
involve partition functions rather than correlation functions of boundary operators. The exponents appearing in (2) 
have already been derived in the limit of large modulus using a radial version of SLE [11,13], and it would be very 
interesting to use these methods to verify the more detailed results of the present paper. 
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